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ABSTRACT 


This  report  describes  the  evaluation  of  the  free 
energy  of  the  nearest  neighbor  Ising  model  through 
tenth  order  in  the  inverse  temperature  variable 
( J/kT ).  The  free  energy  is  expressed  in  terms  of 
diagrams  composed  of  interaction  bonds  and  ver¬ 
tices.  An  explicit  list  of  the  diagrams  with  10  or  fewer 
bonds  is  presented,  and  the  unrestricted  lattice  con¬ 
stants  (free  multiplicities)  of  these  diagrams  are  given 
for  the  cubic  lattices.  We  employ  a  theorem  that 
states  that  a  certain  sum  of  diagrams  is  zero  in  order 
to  obtain  the  density-dependent  weight  factors 
associated  with  the  vertices.  These  factors  are 
polynomials  in  the  density  and  depend  on  the  kind  of 
articulated  vertex  under  consideration.  Explicit 
denumeration  of  the  diagrams  is  given  for  those  ar¬ 
ticulated  vertices  of  the  fourth  and  fifth  kind  that  are 
needed  in  the  tenth  order  expansion.  These  methods 
reproduce  the  known  results  for  the  cubic  lattices. 
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1.  INTRODUCTION 


The  exact  solution  to  the  free  energy  of  the  three- 
dimensional  Ising  model  is  not  known,  but  exact 
series  expansions  have  proven  quite  beneficial  for 
deducing  the  thermodynamic  properties  of  these 
systems.1 * 3*  The  terms  in  these  series  are  evaluated 
through  some  order  in  an  expansion  variable,  and 
diagrammatic  methods  are  a  valuable  aid  for  catalog¬ 
ing  the  various  contributions.  The  steps  involved  in 
obtaining  such  expansions  are:  (1)  list  the 
topologically  distinct  diagrams,  (2)  generate  the 
weight  factors  that  are  associated  with  the  vertices  of 
the  diagram,  and  (3)  calculate  the  appropriate  lattice 


constant  (the  number  of  ways  the  diagram  can  be 
embedded  in  the  lattice).  The  degree  of  difficulty  of 
these  steps,  as  well  as  the  need  to  make  modifications 
appropriate  to  the  particular  model  (e.g.,  particular 
lattice  or  particular  interaction  potential),  depend  on 
the  specific  diagrammatic  expansion.  The  present 
report  uses  an  especially  convenient  expansion  to 
analyze  nearest-neighbor  models  on  cubic  lattices; 
application  to  other  systems  can  be  achieved  by 
calculating  the  unrestricted  lattice  constants  (i.e., 
free  multiplicities)  for  them. 


1  Phase  Transitions  and  Critical  Phenomena  Vol.  3,  C.  Domb 
and  M.  S.  Green  (eds.).  Academic  Press,  Londor  (1974). 

1  M.  F.  Sykes,  D.  S.  Gaum,  P.  D.  Roberts,  and  J.  A.  Wyles, 
"High  Temperature  Series  for  the  Susceptibility  of  the  Ising 
Model;  I.  Two  Dimensional  Lattices,  and  II.  Three  Dimensional 
Lattices,”  J.  Phys.  AS,  6, 640(1972). 

3  M.  F.  Sykes,  D.  L.  Hunter,  D.  S.  McKenzie,  and  B.  R.  Heap, 
“Specific  Heat  of  a  Three  Dimensional  Ising  Ferromagnet 
above  the  Curie  Temperature  II,"  J.  Phys.  A  5, 667  (1972). 


4  P.  H.  E.  Meijer  and  R.  A.  Farrell,  "Pad<  Approximations  and 
the  Critical  Exponents  in  the  Two-and  Three-Dimensional  Ising 
Models,”  Phys.  Rev.  B 12,  3  (1975). 

5  D.  S.  Gaunt  and  G.  A.  Baker,  Jr.,  “Low-temperature  Critical 
Exponents  from  High-temperature  Series:  The  Ising  Model,” 
Phys.  Rev.  B 1, 1184(1970) 

4  G.  Paul  and  H.  E.  Stanley,  "Partial  Test  of  the  University 
Hypothesis:  The  Case  of  Next-Nearest-Neighbor  Interactions,” 
Phys.  Rev.  BS,  3715(1972). 
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2.  FREE  ENERGY 


The  Helmholtz  free  energy,  F,  of  the  Ising  model 
can  be  written’  *  as 

-$F=  -Uplnp  +  0-p){n(\-p))  +  £  £D', 

S//w  D'  ^  j  ^ 

where  0  is  the  usual  Boltzmann  temperature  factor, 
\/{kT)  ;  p  is  the  density  (fraction  of  sites  that  are  oc¬ 
cupied),  and  D'  is  the  set  of  topologically  distinct 
connected  diagrams  composed  of  interaction  bonds 
and  density  dependent  vertices.  Any  number  of 
bonds  can  connect  sites  i  and  j\  the  weight  to  be 
assigned  such  a  bond  is  (0J:j)  if  i  *  j  and  0  if  i  =  j , 
where  Ju  is  the  interaction  energy  between  particles 
on  these  sites.  The  weight  factor  associated  with  a 
vertex  is  a  polynomial  in  the  density  whose  coeffi¬ 
cients  depend  on  the  “kind  of  articulation  vertex” 
that  is  under  consideration.  An  articulation  vertex  is 
of  the  p\i\  kind  (or  order)  if  its  removal  would  cause 
the  diagram  to  separate  into  p  disjoint  parts.  The 
simplest  example  of  an  articulation  vertex  is  the  cen¬ 
tral  vertex  in  a  star  diagram  composed  of  dumbbells 
that  are  joined  at  (or  share)  a  common  vertex  (see 
Fig.  I).  Such  a  vertex  is  of  the  sth  kind  if  there  are  s 
points  to  the  star. 


Figure  1  —Simple  examples  of  diagrams  with  articulated 
vertices:  »>1. 


The  usual  virial  expansion  of  the  free  energy  in¬ 
volves  the  sum  of  more  than  singly  connected  single- 
bonded  diagrams  composed  of  vertices  having  weight 
p  interconnected  by  bonds  btJ  and  -6 u.  Here, 
single-bonded  means  that  no  more  than  one  bond  can 
connect  a  pair  of  vertices;  b,j  =  [txp(0J:J )  -  IJ; 
and  6,j  is  I  if  /  =  j,  and  is  0  if  /  ytj.  Equation  1  is  ob¬ 
tained’*  by:  (1)  Taylor  series  expanding  the  factor 
[cxp{0Ju )  -  1),  which  leads  to  multiple  (0J  ,,)- 
bonds  and  a  symmetry  factor  to  account  for  the  in¬ 
terchangeability  of  the  multiple  bonds  that  connect  a 
pair  of  vertices;  (2)  using  the  well-known  result  for 
the  hard  core  gas  to  evaluate  the  sum  of  diagrams 
that  do  not  contain  (0JU)- bonds;  (3)  demonstrating 
that  there  is  an  exact  cancellation  among  the  set  of 
diagrams  that  contain  b,d -bonds  and  have  one  or 
more  vertices  that  do  not  have  at  least  one  b;/-bond 
emanating  from  them;  and  (4)  introducing  density 
weight  factors  to  accomplish  a  vertex  renormaliza¬ 
tion,  i.e.,  to  sum  over  -6v-bond  configurations. 
The  density  weight  factor  for  an  articulation  vertex 
of  the  pth  kind  is  denoted  by  Af„(  % ,  where  n,  is 
the  number  of  bonds  that  connect  the 'vertex  to  the 
/th  disjoint  part  into  which  the  diagram  would 
separate  if  the  vertex  were  removed.  Step  3  demon¬ 
strates  that  each  vertex,  p ,  has  one  or  more  b,d  bonds 
emanating  from  it,  so  that  the  weight  factor  is  a  poly¬ 
nomial  in  the  density  of  order  of  the  number  of 

v 

bonds  emanating  from  the  circle,  i.e.,  of  order  £  n, . 

i =  I 

The  density  weight  factor  for  articulation  vertices 
of  the  first,  second,  and  third  kind  have  been  ob¬ 
tained*  from  recurrence  relations  among  their  coeffi¬ 
cients.  A  detailed  description  of  the  derivation  of  the 
recurrence  relations  is  given  in  Ref.  7.  The  generating 
functions  that  lead  to  these  relationships  are  a  direct 
consequence  of  a  diagrammatic  theorem’  that  ex¬ 
presses  the  density  weight  factors  for  vertices  of 
higher  kind  in  terms  of  sums  of  products  of  density 
weight  factors  of  lower  kinds.  This  theorem,  or  zero 
sum  rule,  states  that  the  contribution  of  the  diagram 

representing  A f„ldlJ . n>  plus  the  contributions  of  all 

its  associated  diagrams  sum  to  zero.  The  associated 


’  R.  A.  Farrell,  T.  Morita,  and  P.  H.  E.  Meijer,  “Cluster  Expan¬ 
sion  for  the  Ising  Model, ’ ’  J.  Chem.  Phys.  49, 34  9  (1 966). 

*  P.  Esfandiari,  P.  H.  E.  Meijer,  R.  A.  Farrell,  and  S.  Favin, 
"New  Generating  Functions  and  Results  for  the  Density 
Polynomials  of  the  Lattice  Gas,”  Phys.  Rev.  B  24,  1298  (1981), 


also,  “l  rratum:  New  Generating  Functions  and  Results  for  the 
Density  Polynomials  of  the  Lattice  Gas,”  Phys.  Rev,  B  25, 6030 
(1982). 

*  T.  Morita  and  K.  Hiroike,  “A  New  Approach  to  the  Theory  of 
Classical  Fluids,  III  —  General  Treatment  of  Clasical  Systems," 
Prog.  Theor.  Phys.  25,  537  (1961) 
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diagrams  are  constructed  as  follows.  The  v  disjoint  have  used  analogous  classifications  to  construct  the 

parts  into  which  the  diagram  would  separate  if  this  sets  of  associated  diagrams  for  each  of  the  density 

circle  were  removed  are  represented  by  v  labeled  weight  factors  of  the  fourth  and  fifth  kind  needed  to 

black  triangles  on  a  field  (compare  Fig.  2).  Open  calculate  the  free  energy  through  tenth  order  in  in- 

p-gons  (p  =  2  is  a  circle,  p  =  3  a  triangle,  etc.)  are  in-  verse  temperature.  The  use  of  specific  n,’ s  simplifies 

serted  between  the  black  triangles  in  all  possible  com-  the  procedure  by  eliminating  certain  generic 

binations.  Inside  each  p-gon  are  one  or  more  black  diagrams  that  would  occur  for  larger  n,'s.  In  order  to 

dots  that  are  connected  to  the  black  triangles  by  illustrate  this  simplification,  consider  what  would 

bonds  of  weight  1.  The  following  three  rules  must  be  happen  if  we  had  listed  the  associated  diagrams  for 

obeyed:  (1)  each  black  triangle  is  not  articulated,  and  the  v  =  3  case  of  Af2  2  2  .  Because  only  two  bonds 

black  triangle  /  has  exactly  n,  bonds  emanating  from  would  emanate  from  each  star,  neither  the  second 

it,  (2)  each  black  dot  within  a  given  p-gon  is  con-  nor  the  third  types  of  generic  diagrams  in  Fig.  2 

nected  to  p  specific  black  triangles,  and  (3)  the  would  appear  in  the  list.  Analogous  exclusions  occur 

number  of  black  dots  within  a  p-gon  is  only  restricted  for  v  ~  4  and  5.  In  Ref.  7,  the  procedure  was  il- 

by  the  constraint  that  it  obey  these  rules.  lustrated  by  listing  the  diagrams  for  the  case  of 

This  classification  is  trivial  for  v  =  2,  and  for  v  =  Af2  2222 .  In  this  report,  we  list  all  of  the  required 

3  it  generates  the  five  types  of  diagrams  shown  in  diagrams  in  Figs.  3  through  7.  From  these  figures  and 

Fig.  2.  The  evaluation  of  these  generic  diagrams  pro-  the  generating  functions  derived  in  Ref.  7,  one  ob- 

duces  the  generating  functions  given  in  Ref.  7.  We  tains  the  density  weight  factors  listed  in  Table  1. 


i 


Figure  2— The  five  generic  diagrams  generated  by  the  »  =  3  sum  rule. 
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3m2.2M2  3M|i2/2  6M22M^  6M|M 


Figure  3— The  diagrams  that  contribute  to  <W2  22  2 
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"4-  A 


M3  2  2  2^ 


M5M4/2  M3.3M3/6  3M4jM3/4  3M32M4/4 


M3  j 


3M3  ,  ,M2/8  3M,  ,2M3/8  3M3  2M§/2  3M3  2M|/2 


3M22M3Mj/2  3M4M3M 


M3,2M2,2/4  3M32M22/4 


I 


M3  2M|/2  3M3  jM|/2  3M2  2M3M2/2  3M,  2M3Mj/2  M3! 


6M3M^ 


3M4M3M2  3M4M3M2/2 


3Mj  jMjMj 


Figure  4— The  diagrams  that  contribute  to  M3  2>2,2. 
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Figure  5— The  diagrams  that  contribute  to  M4  2.2,2  ■ 
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Table  1  — The  coefficients  W„v„  m  In 
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3.  APPLICATION  TO  THE  NEAREST  NEIGHBOR  MODEL 


In  order  to  determine  the  free  energy  we  con¬ 
structed  two  sets  of  lattice  diagrams.  The  first  set 
consists  of  ail  connected  and  nonarticulated 
diagrams  that  have  10  or  fewer  bonds.  There  are  185 
single-bonded  diagrams  of  this  type;  they  are  listed  in 
Fig.  8,  together  with  their  symmetry  numbers.  The 
free  multiplicities  of  these  diagrams  for  nearest- 
neighbor  models  on  the  cubic  lattices  are  presented  in 
Tables  2  and  3.  They  are  in  accord  with  the  recently 
published  tables  of  Kincaid,  Baker,  and  Fullerton. 10 
The  exclusion  of  triangular,  pentangular,  etc.  figures 
on  the  simple-cubic  and  body-centered-cubic  lattices 
greatly  reduces  the  number  of  diagrams  needed  for 
those  cases. 

The  second  set  of  lattice  diagrams  consists  of  all 
the  articulated  diagrams  with  10  or  fewer  bonds. 
These  are  generated  by  joining  together  two  or  more 
nonarticulated  diagrams  at  one  or  more  common  ver¬ 
tices.  In  constructing  these  diagrams,  one  uses  all  the 
previously  listed  nonarticulated  diagrams  except  that 
the  single-bonded  two-vertex  diagram  (i.e.,  the 
single-bonded  dumbbell)  cannot  occur  as  an  element 
emanating  from  a  vertex;  such  a  diagram  element  has 
to  be  at  least  double-bonded.  This  exclusion  is  a  con¬ 
sequence  of  the  fact  that  the  articulation  vertices  are 
formed  by  collapsing  a  number  of  vertices  (that  were 
connected  by  -  6  ,y -bonds).  Subject  to  this  restric¬ 
tion,  we  Find  that  there  are  245  topologically  distinct 
articulated  diagrams  with  the  minimum  number  of 
bonds  between  their  vertices;  they  are  listed  in  Fig.  9, 
together  with  their  symmetry  numbers.  The  free 
multiplicity  of  an  articulated  diagram  is  simply  the 
product  of  the  free  multiplicities  of  the  nonar¬ 
ticulated  diagram  elements  used  to  generate  it;  it  has 
not  been  listed. 

After  these  lists  were  composed  and  checked,  we 
generated  the  additional  terms  by  multiple-bonding 
all  of  the  previously  listed  diagrams.  Multiple¬ 
bonding  does  not  change  the  number  of  ways  the 
diagram  fits  into  the  lattice  (i.e.,  its  free  multiplicity) 
but  it  does  change  the  density  weight  factors,  the 
order  of  the  temperature  contribution,  and  the  sym¬ 
metry  number  appropriate  to  the  diagram.  Multiple¬ 
bonding  changes  the  symmetry  of  a  diagram  for  two 
reasons.  First,  multiple-bonds  that  connect  a  given 
pair  of  vertices  are  equivalent;  this  is  accounted  for 


ln  J.  M.  Kinkaid,  G.  A.  Baker,  Jr.,  and  L.  W.  Fullerton,  "High- 
Temperature  Series  Expansions  of  the  Continuous-Spin  Ising 
Model,"  Los  Alamos  Scientific  Report  La-Ur-79-1575  (un¬ 
published). 


by  introducing  a  product  of  factorials  in  recognition 
of  their  interchangeability.  The  second  reason  is  that 
the  overall  symmetry  of  the  single-bonded  diagram 
can  be  changed.  For  example,  the  symmetry  number 
of  a  p-gon  is  2 p  because  there  are  p  rotations  and  one 
flip  of  the  diagram,  each  of  which  leads  to  a  diagram 
that  is  topologically  equivalent.  If  one  of  the  pairs  of 
connected  points  in  a  p-gon  is  joined  by  multiple 
bonds  and  the  other  (p-  I)  connections  are  single- 
bonded,  then  the  symmetry  number  is  reduced  to  2 
because  the  rotations  lead  to  topologically  distinct 
diagrams;  only  the  flip  around  the  center  of  the 
multiple  bonded  pair  leads  to  a  topologically 
equivalent  diagram.  Usually  one  would  account  for 
this  change  of  symmetry  by  listing  only  topologically 
distinct  unlabeled  diagrams  and  using  their  modified 
symmetry  numbers.  Alternatively,  one  can  sum  over 
all  possible  multiple  bondings  of  the  single-bonded 
diagram  independently  arid  use  the  overall  symmetry 
number  of  the  single-bonded  diagram  (modified  by 
the  product  of  factorials  mentioned  above).  We  have 
used  this  alternative  because  it  is  ideally  suited  for  a 
computing  machine.  This  procedure  is  based  on  the 
fact  that  independent  multiple-bonding  can  lead  to 
topologically  equivalent  diagrams,  and  the  number 
of  such  diagrams  of  a  given  type  is  simply  the  ratio  of 
the  symmetry  number  of  the  single-bonded  diagram 
to  that  of  the  multiple  bonded  diagram.  A  few  simple 
examples  are  presented  in  Fig.  10. 

We  created  the  computer  compilation  program 
reproduced  and  annotated  in  the  appendix  in  order  to 
accumulate  the  contributions  of  the  aforementioned 
430  lattice  diagrams  and  the  needed  multiple-bonded 
diagrams  that  they  generate.  The  following  informa¬ 
tion  is  used  to  characterize  each  diagram:  an  ar¬ 
bitrary  label,  1  to  430;  the  symmetry  number  of  the 
single-bonded  diagram;  the  unrestricted  lattice  con¬ 
stant  (i.e.,  free  multiplicity)  for  the  lattice  under  con¬ 
sideration;  the  total  number  of  vertices,  NC\  the  total 
number  of  bonds  (actually,  the  number  of  connected 
pairs  of  vertices),  NB,  in  the  single-bonded  diagram; 
and  two  sequences  of  numbers  to  specify  the 
diagram’s  topology  (i.e.,  the  manner  in  which  the 
vertices  are  interconnected).  The  first  sequence  gives 
the  lower  limits  on  the  number  of  bonds  between 
each  of  the  NB  pairs  of  vertices  that  are  connected. 
These  numbers  are  in  general  equal  to  unity,  except 
in  the  case  of  a  two-vertex  element  (dumbbell)  in  an 
articulated  diagram  where  the  number  is  2  (for  the 
reasons  discussed  above).  The  second  sequence  of 
numbers  consists  of  two  subsets.  The  first  subset 


21 


FriSCXOIMO  BttM  MAW-MOT  TILMD 


TMt  JOHN*  HOPKINS  UNMIWTV 

APPLIED  PHYSICS  LABORATORY 

IAUML.  MAOTUMO 


specifies  the  articulation  order  of  each  of  the  NC  ver¬ 
tices  ;  i.e.t  it  is  a  sequence  of  NC  numbers  that  are 
unity  if  the  vertex  is  nonarticulated,  2  if  the  articula¬ 
tion  is  of  the  second  kind,  etc.  The  second  subset 
gives  the  emanation  numbers  of  each  vertex.  If  the 
vertex  is  nonarticulated,  there  is  only  one  set  of 
emanation  numbers  that  specify  the  other  vertices  to 
which  it  is  connected.  If  it  is  an  articulation  vertex  of 
the  pth  kind,  there  are  p  sets  of  emanation  numbers 
needed  because  the  diagram  separates  into  p  elements 
if  the  vertex  is  removed.  In  this  case,  each  set  of 
emanation  numbers  specifies  those  vertices  within  a 
particular  one  of  the  p  elements  that  are  connected  to 
the  vertex  of  interest.  This  scheme,  cumbersome  as  it 
may  seem,  is  a  very  efficient  way  to  transfer  the  in¬ 
formation  into  the  code.  The  data  were  proofread 
several  times  and  were  checked  by  using  them  to 
regenerate  the  diagram. 

The  program  takes  these  data  for  a  single-bonded 
diagram,  generates  the  labels  of  the  appropriate  den¬ 
sity  weight  factors,  multiplies  these  NC  polynomials 
together,  and  multiplies  the  result  by  the  ratio  of  the 
free  multiplicity  divided  by  the  symmetry  number. 
The  final  result  is  then  accumulated  in  the  appropri¬ 


ate  Cnm,  where 

-  0F  =  -  L[p< np  +  (1  -  p)  f n(l  -  p)J 
10  2* 

+  L  L  C„,mpm(J/kD”  +  OHJ/kT )"J. 

«-l  m-l  (2) 

Next,  the  program  automatically  multiple-bonds  the 
diagram  and  repeats  the  procedure.  The  free  energy 
series  for  the  three  cubic  lattices  are  given  in  Tables  4 
through  6.  These  coefficients  were  compared  with  the 
results  available  in  the  literature.  In  particular,  for 
the  simple-cubic  and  body-centered-cubic  lattices, 
the  low  temperature  series  of  Sykes,  Essam,  and 
Gaunt11  were  used  to  determine  the  high  temperature 
series  by  means  of  an  inversion  of  the  activity  series 
into  a  density  series.  For  the  face-centered-cubic  lat¬ 
tice,  McKenzie’s  recent  results'2  for  the  low  tempera¬ 
ture  series  were  inverted.  Exact  agreement  is  found  in 
all  cases.  Thus,  we  are  confident  that  the  symmetry 
numbers  and  density  weight  factors  are  correct.  The 
evaluation  of  Eq.  1  for  more  general  models  (e.g., 
longer  ranged  interactions,  multiple  sublattice  long 
range  order  parameters,  etc.)  are  under  way. 


M.  F.  Sykes,  J.  W.  Essam,  and  D.  S.  Caunt,  "Derivation  of  13  S.  McKenzie,  "Extended  High  Temperature  Low  Field  Expan- 
Low-Temperature  Expansions  for  the  Ising  Model  of  a  Fer  ro-  sions  for  the  I  sing  Model,"  Can.  J.  phys.  51. 1239(1979). 
magnet  and  an  Antiferromagnet,  ”  J.  Math.  Phys.  4, 283  ( 1965). 
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Table  2— The  free  multiplicities  for  the  simple-cubic  and  body-centered-cubic  lattices. 


FM 

FM 

FM 

FM 

JL 

sc 

bcc 

L 

sc 

bcc 

1 

6 

8 

105 

18612 

140608 

3 

90 

216 

106 

18612 

140608 

7 

318 

1000 

109 

20250 

157464 

9 

1860 

8000 

115 

8166 

39304 

14 

1350 

5832 

128 

20250 

157464 

22 

1494 

5832 

129 

18324 

140608 

24 

960 

4096 

131 

16134 

125000 

30 

5712 

32768 

132 

20250 

157464 

32 

44730 

343000 

135 

129330 

1331000 

34 

1172556 

16003008 

137 

108600 

1124864 

36 

27900 

216000 

186 

36 

64 

59 

642 

2744 

187 

216 

512 

64 

3804 

21952 

188 

216 

512 

67 

4770 

27000 

189 

1296 

4096 

69 

22632 

175616 

190 

1296 

4096 

65 

4044 

21952 

191 

1296 

4096 

86 

3102 

17576 

192 

7776 

32768 

103 

24216 

175616 

193 

7776 

32768 

phsckdiio  not  kjjb«not  nue 


25 


THE  JOHNS  HOPKINS  UNIVERSffY 

APPUED  PHYSICS  LABORATORY 

LAUREL  Mary  LAND 


Tablv  3— The  free  multiplicities  for  the  face  centered-cubic  lattice. 
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Figure  10— Examples  of  the  two  alternate  ways  ot  multiple-bonding  a  diagram.  (As  discussed  in  the  text,  the  first  method 
uses  only  topologically  distinct  diagrams  and  their  associated  symmetry  numbers.  The  alternative  procedure,  which  was 
used,  does  the  multiple-bonding  independently  and  uses  only  the  symmetry  number  of  the  single-bonded  diagram.  The  advan¬ 
tage  of  this  alternative  is  that  it  is  trivial  to  implement  on  the  computer.) 
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Table  4— The  coefficients  of  />*  [i^J)n  in  the  Tayior-series  expansion  of  the  Helmholtz  free  energy  of  the  nearest-neighbor  lat¬ 
tice  gas  on  the  face-centered-cubic  lattice.  The  entry  C<  n.A )  =  //./mis  to  be  read  as//  +  [/m/(n'2)],  e  g.,  C(9,8)  =  21591:498960 
=  215291  +  498960/(9!  2)  =  215291.6875. 


k*/rv»l  23456  7  8  9  10 


2 

6:0 

3:0 

1  :0 

0:  12 

0:  12 

0:12 

0:  12 

0:12 

0:  12 

0:1? 

3 

0:0 

-6:0 

2:0 

8:24 

8:  120 

5:696 

2: 7224 

1 : 9048 

0 : 284280 

0:663416 

4 

0:0 

3:0 

-11:0 

-34:  12 

-10:160 

49:372 

96:586 

104:66036 

84:613740 

56:615252 

5 

0:0 

0:0 

12:0 

57:0 

-75:0 

-374:720 

-  394:7056 

158:17136 

1 042 : 344860 

1765:2409640 

6 

0:0 

0:0 

-4:0 

-54:0 

?65 : 96 

1021:0 

-149:4366 

-4336:14112 

-7039:495964 

-1164: 8724000 

7 

0:0 

0:0 

0:0 

30:0 

-468:0 

-1495:0 

4542:0 

16482:30240 

-6581 : 181440 

-114154:5320600 

6 

0:0 

0:0 

0:0 

-7:24 

423:0 

1225:1060 

-14813: 7920 

-21321:7560 

215291 : 498*60 

916676:3670360 

9 

0:0 

0:0 

0:0 

0:0 

-204:0 

-460:0 

26654:0 

-33913:40320 

-  1  119307:544320 

-4000946 : 3626600 

10 

0:0 

0:0 

o.o 

0:0 

40:192 

8:0 

-31784:0 

19324  7 . 0 

3419749:84672 

11976410:4596460 

1 1 

0:0 

0:0 

0:0 

0:0 

0:0 

46:0 

25452:0 

-393420.0 

-7211599:0 

-26566250: 1451520 

12 

0:0 

0:0 

0:0 

0:0 

0:0 

-6:0 

-13446:0 

4 768900 

11052791: 362880 

45106674:5564160 

13 

0-0 

0:r 

0:0 

0:0 

0:0 

0:0 

4248:0 

-371007:0 

-12494314:0 

-59478333:0 

14 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

-606:6640 

182406:0 

10350458:0 

61 194742:0 

15 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

-51762:0 

-6122232:0 

-48994063:4354560 

16 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

6470:20160 

2450285:0 

30169646: 1088640 

17 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

-594708:0 

-13920714:0 

16 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

66078:483840 

4565694:0 

19 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

-953352:0 

20 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

0:0 

95335:1451520 

Table  5— The  coefficients  of  jj*  (/i/)n  in  the  Tayior-series  expansion  of  the  Helmholtz  free  energy  of  the  nearest-neighbor  lat¬ 
tice  gas  on  the  body-centered  cubic  lattice 
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Tabla  8— The  coefficients  of  UiJ)n  in  the  Tayior-series  expansion  of  the  Helmholtz  free  energy  of  the  nearest-neighbor  lat¬ 
tice  gas  on  the  simple  cubic  lattice. 
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APPENDIX 


THE  COMPILATION  PROGRAM 


The  notations  used  in  this  computer  program  are: 


1 


Diagram  label. 


quent  labels  are  \jao(J)  -  1J  oui 
of  sync.  (This  concept  is  il¬ 
lustrated  below.) 


S 

W 


Symmetry  number. 
Free  multiplicity. 


NP  — (5,NNN)  Specifies  the  bonds  that  ema¬ 
nate  from  the  vertex  labeled  by 
NNN.  (Again,  see  below.) 


NB 


Number  of  connected  vertex  11T 
pairs. 


L  —  L(IO);  L(l),  The  lower  limit  on  the  number 
L(2).  ...,L(NB)  of  bonds  between,  respectively,  IjT 

the  first,  second . and  NBih 

vertex  pair. 


Specifies  the  total  number  of 
bonds  that  can  be  added  to  a 
single-bonded  diagram. 

Specifies  the  upper  limit  on  the 
number  of  bonds  between  the 
yth  vertex  pair. 


INSLUM  The  number  of  bonds  in  the 

single-bonded  diagram. 

NC  The  number  of  vertices  (circles). 

JAO  —  JAO(NC);  The  articulation  orders  of,  re- 

JAOfl),  JAO(2)  spectively,  the  first,  second . 

...,JAO(NC)  and  /VCth  vertex. 


JAOSUM  The  number  of  indices  that 

must  be  specified  in  the  product 
of  polynomials  (Af  s). 

NNN  Specifies  the  vertex  label  unless 

one  (or  more)  of  the  vertices 
(say  the  yth)  is  articulated,  in 
which  case  jao(j)  of  the  NNN ’s 
are  used  for  it  and  all  subse- 


1NL  —  L(N). 

Figures  A-l  and  A-2  illustrate  these  conventions 
for  two  representative  articulated  diagrams. 

Lines  2  through  59  of  the  main  program  read  in  the 
data,  identify  the  topology  of  the  single-bonded  dia¬ 
gram  and  prepare  the  do-loops  to  be  used  in  multiple 

bonding.  The  integers  y, ,  y'2 . yN8  generated  in  lines 

61  through  125  specify,  respectively,  the  number  of 
bonds  between  the  first,  the  second,...,  and  the  NBih 
vertex  pairs.  Lines  126  through  168  obtain  the  in¬ 
dices,  N(III),  for  each  of  the  circles.  ///;  identify  the 
corresponding  density  weight  factors;  and  multiplies 
them  together.  The  multiplication  of  these  poly¬ 
nomials  is  accomplished  in  a  subroutine,  POLYMT, 
which  is  also  listed.  Steps  169  through  214  accumu¬ 
late  the  contributions  to  the  coefficients  C„,„  of 
Eq.2. 

The  following  pages  list  the  compilation  program. 


FttECEDUNQ  PAG*.  HUNK -NOT  HLMH> 


33 


no  n  non  non 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL  MARYLAND 


EVALUATES  ARTICULATION  ISING  DIAGRAMS.  S  FAVIN  AUG  1982 

MAIN 

10 

IMPLICIT  REAL*16 ( A-H , 0-2 ) 

MAIN 

20 

DIMENSION  P(ll,3,3,4,6,ll),  L(10),  1(10),  Sl(21) 

MAIN 

30 

DIMENSION  NP(5 , 10) ,  JAO(IO),  N(5) 

MAIN 

40 

DIMENSION  SSI (21, 11) 

MAIN 

50 

DIMENSION  ZZ(26136) 

MAIN 

60 

COMMON  DEBUG 

MAIN 

70 

LOGICAL  DEBUG 

MAIN 

80 

EQUIVALENCE  (P(l , 1 , 1 , 1 , 1 , 1 ) ,  ZZ(1)) 

MAIN 

90 

EQUIVALENCE  (N(1),N1),  (N(2),N2),  (N(3),N3),  (N(4),N4), 

(N(5) ,N5) 

MAIN 

100 

EQUIVALENCE  (L(1),I1L),  (L(2),I2L),  (L(3),I3L),  (L(4),I4L), 

MAIN 

110 

1  (L(5) , I5L) ,  (L(6) , I6L) ,  (L(7),I7L),  (L(8),I8L),  (L(9) 

,I9L), 

MAIN 

120 

2  (L(10) , I10L) 

MAIN 

130 

DATA  SSI  /231*0.0Q0/ 

MAIN 

140 

DATA  L/10*0/ ,  1/10*0/ 

MAIN 

150 

MAIN 

160 

DEBUG  =  .TRUE. 

MAIN 

170 

DEBUG  =  .FALSE. 

MAIN 

180 

DEBUG  =  .FALSE. 

MAIN 

190 

DO  3  K=l, 26136 

MAIN 

200 

3  ZZ(K)  =  0.0D0 

MAIN 

210 

CALL  READIN(  P  ) 

MAIN 

220 

MAIN 

230 

28-ELEMENTAL  (SINGLE-BONDED)  DIAGRAMS 

MAIN 

240 

PRINT  17 

MAIN 

250 

MAIN 

260 

DO  100  J=1 ,430 

MAIN 

270 

IF (DEBUG)  PRINT  17 

MAIN 

280 

17  FORMAT ( 1 1 1 ) 

MAIN 

290 

READ  (5,1 ,END=101)  S,W,NB,NC,L 

MAIN 

300 

1  FORMAT(2F10. 0,1212) 

MAIN 

310 

PRINT  2,  J,  S,  W,  NB,  NC,  L 

MAIN 

320 

2  FORMAT ( 14, 2F16.1, 216, 3X, 1013) 

MAIN 

330 

IF(S  .EQ.  0.0D0)  GO  TO  789 

MAIN 

340 

MAIN 

350 

INLSUM  =  0 

MAIN 

360 

DO  15  M=1 ,NB 

MAIN 

370 

INLSUM  =  INLSUM  +  L(M) 

MAIN 

380 

15  CONTINUE 

MAIN 

390 

MOW  CHECK  ON  NC  1"S  IN  JAO 

MAIN 

400 

MAIN 

410 

READ  5, JAO 

MAIN 

420 

5  FORMAT(10I2) 

MAIN 

430 

IF (DEBUG)  PRINT  6,  ( JAO(M) ,M=1 ,NC) 

MAIN 

440 

PRINT  6,  ( JAO(M) ,M=1 ,NC) 

MAIN 

450 
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6 

FORMAT ( 1 0  JAOs  ',1014) 

MAIN 

460 

J AO SUM  =  0 

MAIN 

470 

DO  16  M=1 ,NC 

MAIN 

480 

JAOSUM  =  JAOSUM  +  JAO(M) 

MAIN 

490 

16 

CONTINUE 

MAIN 

500 

c 

IF (DEBUG)  PRINT  7,  JAOSUM 

MAIN 

510 

PRINT  7,  JAOSUM 

MAIN 

520 

7 

FORMAT ( 1  JAOSUM  IS  ',15) 

MAIN 

530 

c 

MAIN 

540 

c 

ASSUME  JAOSUM  >  0 

MAIN 

550 

DO  18  NNN=1, JAOSUM 

MAIN 

560 

READ  5,  (NP(JJ,NNN) , JJ=1 , 5) 

MAIN 

570 

c 

IF (DEBUG)  PRINT  8,  NNN,  (NP( JJ ,NNN) , JJ=1 , 5) 

MAIN 

580 

PRINT  8,  NNN,  (NP( JJ ,NNN) , JJ=1 , 5) 

MAIN 

590 

18 

CONTINUE 

MAIN 

600 

8 

FORMAT(I5 , 513) 

MAIN 

610 

c 

MAIN 

620 

c 

MAIN 

630 

c 

MAIN 

640 

c 

MAIN 

650 

IIT  =  10  -  INLSUM 

MAIN 

660 

IIT  =  IIT  +  I1L 

MAIN 

670 

I2T  =  IIT  +  *2L 

MAIN 

660 

I3T  =  IIT  +  I3L 

MAIN 

690 

I4T  *  IIT  +  I4L 

MAIN 

700 

I5T  *  IIT  +  I5L 

MAIN 

710 

I6T  =  IIT  +  I6L 

MAIN 

720 

I7T  =  IIT  +  I7L 

MAIN 

730 

I8T  =  IIT  +  I8L 

MAIN 

740 

I9T  =  IIT  +  I9L 

MAIN 

750 

HOT  =  IIT  +  I10L 

MAIN 

760 

10 

FORMAT ( '0  IIT= 1 , 15 , 5X, 'IIT, . . . I10T= ’ , 1015) 

MAIN 

770 

IF(DEBUG) 

MAIN 

780 

*PRINT  10,  IIT,  I1T , I2T , I3T, I4T , I5T , I6T , I7T , I8T , I9T , HOT 

MAIN 

790 

c 

MAIN 

800 

DO  90  I1=I1L,  IIT 

MAIN 

810 

I(1)=I1 

MAIN 

820 

IF(NB  .EQ.  1)  GO  TO  999 

MAIN 

830 

DO  88  I2=I2L,  I2T 

MAIN 

840 

J2  =  12  +  11 

MAIN 

850 

IF(J2  .GT.  10)  GO  TO  90 

MAIN 

860 

1(2)  =  12 

MAIN 

870 

IF(NB  .EQ.  2)  GO  TO  999 

MAIN 

680 

DO  86  I3*I3L,  I3T 

MAIN 

890 

J3  *  J2  +  13 

MAIN 

900 
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IF(J3  .GT.  10)  GO  TO  88 
1(3)  =  13 

IF(NB  .EQ.  3)  GO  TO  999 
DO  84  I4=I4L ,  I4T 
J4  =  J3  +  14 
IF(J4  .GT.  10)  GO  TO  86 
1(4)  =  14 

IF (NB  .EQ.  4)  GO  TO  999 
DO  82  I5=I5L,  I5T 
J5  =  J4  +  15 
IF(J5  .GT.  10)  GO  TO  84 
1(5)  =  15 

IF (NB  .EQ.  5)  GO  TO  999 
DO  80  I6=I6L,  I6T 
J6  =  J5  +  16 
IF(J6  .GT.  10)  GO  TO  82 
1(6)  =  16 

IF (NB  .EQ.  6)  GO  TO  999 
DO  78  I7=I7L,  I7T 
J7  =  J6  +  17 
IF(J7  .GT.  10)  GO  TO  80 
1(7)  =  17 

IF (NB  .EQ.  7)  GO  TO  999 
DO  76  I8=I8L, I8T 
J8=J7+I8 

IF(J8  .GT.  10)  GO  TO  78 
1(8)  =  18 

IF(NB  .EQ.  8)  GO  TO  999 
DO  74  I9=I9L, I9T 
J9  =  J8  +  19 
IF(J9  .GT.  10)  GO  TO  76 
1(9)  =  19 

IF(NB  .EQ.  9)  GO  TO  999 
DO  72  I10=I10L ,  HOT 
J10  =  J9  +  110 
IF(J10  .GT.  10)  GO  TO  74 
1(10)  =  110 
C 
C 
C 
C 
C 

999  F  =  S 

IF (DEBUG) 

*PRINT  12,  II,  12,  13,  14,  15,  16,  17 


MAIN  910 
MAIN  920 
MAIN  930 
MAIN  940 
MAIN  950 
MAIN  960 
MAIN  970 
MAIN  980 
MAIN  990 
MAIN1000 
MAIN1010 
MAIN1020 
MAIN1030 
MAIN1040 
MAIN1050 
MAIN1060 
MAIN1070 
MAIN1080 
MAIN1090 
MAIN1100 
MAIN1110 
MAIN1120 
MAIN1130 
MAIN1140 
MAIN1I50 
MAIN1160 
MAIN1170 
MAIN1180 
MAIN1190 
MAIN1200 
MAIN1210 
MAIN1220 
MAIN1230 
MAIN1240 
MAIN1250 
MAIN1260 
MAIN1270 
MAIN1280 
MAIN1290 
MAIN1300 
MAIN1310 
MAIN1320 
MAIN1330 
MAIN1340 

18,  19,  110  MAIN1350 
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12  F0RMAT{'0  11 . 110=  ',1015) 

MAIN1360 

DO  20  K=2 ,21 

KAIN1370 

20 

Si (K)  =  0.0D0 

MAIN1380 

Sl(l)  =  1.0D0 

MAIN1390 

NS1  =  1 

MAIN1400 

NNN=0 

MAIN1410 

DO  30  11=1, NC 

MAIN1420 

I AO  =  JAO(II) 

MAIN1430 

DO  29  111=1,5 

MAIN1440 

29 

N(III)=0 

MAIN1450 

DO  27  III  =  1,  I AO 

MAIN1460 

NNN=NNN+1 

MAIN1470 

DO  26  JJ=1,5 

MAIN1480 

KK=NP(JJ,NNN) 

MAIN1490 

IF(KK  . LE.0)  GO  TO  27 

MAIN1500 

N(III)  =  N(III)  +  I (KK) 

MAIN1510 

26 

CONTINUE 

MAIN1520 

27 

CONTINUE 

MAIN1530 

MAIN1540 

NOW  ORDER  THE  N(I)  AND  MULT  SI  *  P(N,  Nl+1,  N2+1 ,  N3+1) 

MAIN1550 

MAIN1560 

MAIN 1570 

DO  35  11=1,4 

MAIN1580 

IS=IN+1 

MAIN1590 

DO  35  JN=IS , 5 

MAIN1600 

IF(N(IN)  .LE.  N(JN))  GO  TO  35 

MAIN1610 

NSAVE  =  N( IN) 

MAIN1620 

N( IN)  =  N(JN) 

MAIN1630 

N(JN)  =  NSAVE 

MAIN1640 

35 

CONTINUE 

MAIN1650 

IF (DEBUG)  PRINT  11,  N 

MAIN166 0 

11 

FORMAT ( 1  SORT  N  IS  ' ,515) 

MAIN1670 

Nil  =  Nl+1 

MAIN1680 

N21=  N2+1 

MAIN1690 

N31=  N3+1 

MAIN1700 

N41=  N4+1 

MAIN1710 

N51=  N5+1 

MAIN1720 

MAIN1730 

NN1  =  N1  +  N2  +  N3  +  1  +  N4  +  N5 

MAIN1740 

CALL  P0LYMT(P(1,N11,  N21,  N31,  N41,  N51),  NN1 ,  SI,  NS1 ,  SI) 

MAIN1750 

NS1  «  NS1+NN1-1 

MAIN1760 

30 

CONTINUE 

MAIN1770 

MAIN1780 

MAIN1790 

MAIN180Q 
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KK=0 

MAIN1810 

MAIN1820 

MAIN1830 

F=S 

MAIN1840 

DO  40  K=1 ,NB 

MAIN1850 

FACT  =  1.0D0 

MAIN1860 

IFK  =  I (K) 

MAIN1870 

DO  28  IFF=1 , IFK 

MAIN1880 

28 

FACT  =  FACT  *  IFF 

MAIN1890 

F  =  F  *  FACT 

MAIN1900 

IF (DEBUG) 

MAIN1910 

*PRINT  9,  NB ,NC , I ,  F 

MAIN1920 

9 

FORMAT ( 1  FACT: 1 ,216, 1013, 1PG20. 12  ) 

MAIN1930 

40 

KK  =  KK  +  I (K) 

MAIN1940 

IF(KK  .GT.  0)  GO  TO  49 

MAIN1950 

PRINT  13,  I 

MAIN1960 

13 

FORMAT( 1 l****ERROR  KK=0  AND  I=',10I5) 

MAIN1970 

STOP 

MAIN1980 

49 

CONTINUE 

MAIN1990 

MAIN2000 

KZ=1 

MAIN2010 

DO  50  K=1 , 21 

MAIN2020 

SS1(K,KK)  =  SSI (K , KK)  +  W*S1(K)/F 

MAIN2030 

IF (SSI (K , KK)  .NE.  0.0D0)  KZ=K 

MAIN2040 

50 

CONTINUE 

MAIN2050 

IF(J  .GT.  3)  GO  TO  60 

MAIN2060 

MAIN2070 

4 

FORMAT ( 1  KK= ' , 15 , 1P5G20 . 12 , / , (1P6G20 . 12 ) ) 

MAIN2080 

MAIN2090 

IF(DEBUG)  PRINT  4,  KK ,  (SSI (K , KK) , K=1 . KZ) 

MAIN2100 

PRINT  4,  KK,  (SSI (K,KK) ,K=1 ,KZ) 

MAIN21 10 

4 

FORMAT ( 1  KK=' ,15,/, (1P5G26.16)) 

MAIN2120 

60 

CONTINUE 

MAIN2130 

GO  TO  (90,  88,  86,  84,  82,  80,  78,  76,  74,  72),  NB 

MAIN2140 

72 

CONTINUE 

MAIN2150 

MAIN2160 

74 

CONTINUE 

MAIN2170 

76 

CONTINUE 

MAIN2180 

78 

CONTINUE 

MAIN2190 

80 

CONTINUE 

MAIN2200 

82 

CONTINUE 

MAIN2210 

84 

CONTINUE 

MAIN2220 

86 

CONTINUE 

MAIN2230 

88 

CONTINUE 

MAIN2240 

90 

CONTINUE 

MAIN2250 
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100 

CONTINUE 

MAIN2260 

101 

CONTINUE 

MAIN2270 

DO  954  KK=1 ,11 

MAIN2280 

PRINT  4,  KK,  (SS1(K,KK) ,K=1,21) 

MAIN2290 

PUNCH  25,(SS1(K,KK) ,K=2,21) 

MAIN2300 

25 

FORMAT (1P3G26. 16) 

MAIN2310 

954 

CONTINUE 

MAIN2320 

789 

CONTINUE 

MAIN2330 

STOP 

MAIN2340 

END 

MAIN2350 

SUBROUTINE  READIN  (  P  ) 

IMPLICIT  REAL*16 ( A-H , 0-2) 

DIMENSION  P( 11 , 3 , 3 ,4 , 6 , 11 ) ,  IP(5) 

EQUIVALENCE  (IP(l).N),  (IP(2),L),  (IP(3),K),  (IP(4),J), 
PRINT  1 

1  FORMAT < 1 1  READIN  P(M,N,L,K,J,I) ,M=2 , 11 ) ' , //) 

DO  100  IR=1 , 100 

READ  (9,4)  IP 

4  FORMAT (51 3) 

IF(N  .EQ.  0)  GO  TO  90 

READ (9,5)  (P(M,N,L,K,J,I) ,M=2 ,11) 

PRINT  6,  N,  L,  K , J , I ,  (P(M,N,L ,K , J, I) , M=2 , 11 ) 

6  FORMAT ('O' , 515 , 5F13 . 2 , / , (25X , 5F13 . 2) ) 

100  CONTINUE 
IR=101 

90  CONTINUE 
IR=IR-1 
PRINT  3,  IR 

3  FORMAT ( 1 0  ',15,'  CARDS  WERE  READIN  INTO  POLY') 

5  FORMAT (10F8.0) 

PRINT  2 

2  FORMAT ( ' l 1 ) 

RETURN 

END 


READ  10 
READ  20 
READ  30 
READ  40 
(IP(5) , I)  READ  50 
READ  60 
READ  70 
READ  80 
READ  90 
READ  100 
READ  110 
READ  120 
READ  130 
READ  140 
READ  150 
READ  160 
READ  170 
READ  180 
READ  190 
READ  200 
READ  210 
READ  220 
READ  230 
READ  240 
READ  250 
READ  260 
READ  270 
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SUBROUTINE  POLYMT(P,  NP,  S,  NS,  A) 

IMPLICIT  REAL*16(A-H,0-Z) 

DIMENSION  P(NP) ,  S(NS),  A(21),  B(21) 

COMMON  DEBUG 
LOGICAL  DEBUG 
C 

C  DEBUG  =  .TRUE. 

C  DEBUG  =  .FALSE. 

IF (DEBUG)  PRINT  1,  NP,  (P( I ) , 1  =  1 , 11 ) 

IF(DEBUG)  PRINT  2,  NS,  (S ( I ) , 1=1 , 21 ) 

1  FORMATC  INPUT  P  POLY  NP=',I5,  // , (1P5G20. 12)  ) 

2  FORMAT ( 1  INPUT  S  POLY  NS=',I5,  // , ( 1P5G20 . 12)  ) 

3  FORMAT ( '  OUTPUT  THE  B  POLY ‘ , // , < 1P5G20 . 12)  ) 

C 

DO  10  1=1,21 
10  B(I)  =  0.0D0 
C  ASSUME  NP  .GE.  2 

IF(NS  .GT.  1)  GO  TO  20 
C  SPECIAL  FOR  NS=1 

DO  15  L=1,NP 
15  B(L)  =  S(l)  *  P(L) 

GENERAL  RETURN 

99  CONTINUE 

IF(DEBUG)  PRINT  3,  B 
DO  100  L=1 , 21 
A(L)  =  B(L) 

100  CONTINUE 
RETURN 


20  CONTINUE 

IF(NP  .GT.  NS)  GO  TO  40 
C  FOR  NP  .LE.  NS 
NP1  =  NP-1 
DO  30  L=1,NP1 
LP1  =  L+l 
DO  30  1=1, L 

B(L)  =  B(L)  +  P(I)  *  S(LPl-I) 
30  CONTINUE 
C 

DO  33  L=NP ,NS 
LP1  =  L+l 
DO  33  1=1, NP 


POLY  10 
POLY  20 
POLY  30 
POLY  40 
POLY  50 
POLY  60 
POLY  70 
POLY  80 
POLY  90 
POLY  100 
POLY  110 
POLY  120 
POLY  130 
POLY  140 
POLY  150 
POLY  160 
POLY  170 
POLY  180 
POLY  190 
POLY  200 
POLY  210 
POLY  220 
POLY  230 
POLY  240 
POLY  250 
POLY  260 
POLY  270 
POLY  280 
POLY  290 
POLY  300 
POLY  310 
POLY  320 
POLY  330 
POLY  340 
POLY  350 
POLY  360 
POLY  370 
POLY  380 
POLY  390 
POLY  400 
POLY  410 
POLY  420 
POLY  430 
POLY  440 
POLY  450 
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B(L)  =  B(L)  +  P(D  *  S(LPl-I) 

33  CONTINUE 
C 

NPNS  =  NP+  NS 
NS.l  =  NS+1 
N2  =  NPNS  -  1 
DO  36  L=NS1,  N2 
LNS  =  L  -  NS 
II  =  NPNS  -  L 
DO  36  1=1,11 

B(L)  =  B(L)  +  P(LNS+I )  *  S(NS1  -  I) 
36  CONTINUE 
GO  TO  99 
C 


NP  .GT.  NS 
40  CONTINUE 
NS1  =  NS-1 
DO  50  L=1 ,NS1 
LPl  =  L+l 
DO  50  1=1, L 

B(L)  =  B(L)  +  S(I)  *  P(LP1  -  I) 

50  CONTINUE 

DO  53  L=NS ,  NP 
LPl  =  L+l 
DO  53  1=1, NS 

B(L)  =  B(L)  +  S(I)  *  P(LP1  -  I) 

53  CONTINUE 
C 

NP1  =  NP+1 
NPNS  =  NS  +  NP 
N2  =  NPNS  -  1 
DO  56  L=NPl ,  N2 
LNS  =  NPNS  -  L 
NPL  =  L-NP 
DO  56  1=1, LNS 

B(L)  =  B(L)  +  S(NPL+I)  *  P(NP1  -  I) 
56  CONTINUE 
GO  TO  99 
END 


POLY  460 
POLY  470 
POLY  480 
POLY  490 
POLY  500 
POLY  510 
POLY  520 
POLY  530 
POLY  540 
POLY  550 
POLY  560 
POLY  570 
POLY  580 
POLY  590 
POLY  600 
POLY  610 
POLY  620 
POLY  630 
POLY  640 
POLY  650 
POLY  660 
POLY  670 
POLY  680 
POLY  690 
POLY  700 
POLY  710 
POLY  720 
POLY  730 
POLY  740 
POLY  75 0 
POLY  760 
POLY  770 
POLY  780 
POLY  790 
POLY  800 
POLY  810 
POLY  820 
POLY  830 
POLY  840 
POLY  850 
POLY  860 


41 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

Laurel  Maryland 


S  •  2 

16480  lor  f.c.c. 

*  -  v  540  for  s  e¬ 

ll  728  for  lee 

NB  ■  6 

l  -  12.  1.  t.  I  1.  0.  0.  0  0.  0) 

INSLUM  »  6 

NC  -  5 

JAO  •  il  2.  1  1.  1.  0.  0.  0.  0.  0) 
JAOSUM  *  6 


Ip)  lab*} s  the  vert  -s 
i  labels  the  bond 


Circle 

NNN 

P.'P 

m 

» 

o 

o 

o 

o 

(2) 

2 

o 

o 

o 

o 

3 

2.  b.  0 

0  0 

13) 

4 

2.  3.  0. 

o 

o 

14) 

5 

3  4  0. 

0  0 

(5) 

6 

4.  5.  0. 

o 

o 

HT  - 

4 

1 1 T  = 

6 

I2T  =■ 

5 

I3T  * 

5 

I4T  * 

5 

I5T  *= 

6 

Figure  A-1— An  illustration  of  the  notations  used  in  the  mam  program 
(The  diagram's  symmetry  number  does  not  account  for  the  inter 
changeability  of  multiple  bonds  between  a  connected  pair  of  vertices  ) 
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Figure  A-2— An  example  of  the  notation  used  in  the  mam  program 
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